Neutrino mixing in matter 
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Three-neutrino mixing in matter is studied through a set of evolution equations which are based on 
a rephasing invariant parametrization. Making use of the known properties of measured neutrino 
parameters, analytic, approximate solutions are obtained. Their accuracy is confirmed by compar- 
ison with numerical integration of these equations. The results, when expressed in the elements 
squared of the mixing matrix, exhibit striking patterns as the matter density varies. 
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It is well-established that neutrino mixing is modified 
by the presence of matter [l[ . Their effect has been used 
in the analyses of solar neutrinos, and is expected to im- 
pact those of the supernova neutrinos, when and if they 
become available. Closer to home, there is a plethora of 
long baseline experiments either in operation or in the 
planning stage. For these studies, it is essential to in- 
clude the matter effects in order to understand neutrino 
mixing at the fundamental level. 

In the literature, effort has been devoted to solving 
problems along this line However, the process in- 
volves the complication of the cubic eigenvalue problems, 
and the results are usually far from transparent for a clear 
extraction of the physical implications. 

In this work we derive a set of evolution equations 
for the neutrino parameters, using a rephasing invari- 
ant parametrization which was developed for three- flavor 
quark mixing. The same formalism can be used in the 
neutrino sector, as long as it is used for lepton num- 
ber conserving processes, such as in neutrino oscillation, 
which will be studied here. It will be shown that the 
coupled equations have simple, analytic, solutions which, 
when compared to the complete numerical solutions, are 
quite accurate. 

For the neutrino mixing (PMNS) matrix (V), we adopt 
the parametrization introduced earlier Q. Briefly, with- 
out loss of generality, one can demand detV = -1-1. There 
are then a set of rephasing invariants 



r 



ijk 



VliV2jV3k = Rijk ~ iJ, 



(1) 



where their common imaginary part can be identified 
with the Jarlskog invariant J Q. Their real parts were 
defined as 

(i?123,^231,^312;-Rl32,^213,^32l) = {xi, X2, X3; yi, y2, Vs)- 

(2) 
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These variables are bounded by ±1: —1 < {xi,yj) < +1, 
with yj < Xi for any (?, j). They satisfy two constraints 

detV = {xi +X2+ X3) - {yi + y2 + y-i) = 1, (3) 
{xiX2 + X2X3 + X3X1) - {yiy2 + y2y3 + VsVi) = 0. (4) 



In addition, it is found that 

j2 = X1X2X3 - yiy2y3- 
The {x,y) parameters are related to P by 



(5) 



xi -yi X2- y2 X3 - y3 
W = [l^^jf] = I x3-y2 xi~y3 X2~yi I . (6) 
X2 -ys X3- yi xi - y2 

One can readily obtain the parameters {x, y) from W by 
computing its cofactors, which form the matrix w with 
w'^W = {detW)I, and is given by 



xi +yi X2 + y2 X3 + y3 
w = \ X3 + y2 xi+y3 X2 + yi 
X2 + y3 X3 + yi xi + y2 



(7) 



Physical measurables can always be expressed in terms 
of (x, y). For instance, the Vfj_ — >■ i/g transition probability 
is given by 
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where Di = neutrino mass squared, L is the length of 
baseline, E is the neutrino energy, and 



X2X3 + xiy2 - yiy2 - y^yz 
-X1X3 - X2X3 + X3yi + ^22/3 
X1X3 + X2y3 - ViVs - y2y3 



^21 
^31 
-'^32 



(9) 



Experimentally, the PMNS matrix in vacuum is well- 
approximated by 



6, 

l+2e^+g 
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-2r, 
-/3 + 77 



(10) 

with (e, ri,f3,£,) ^ 1. reduces to the tri-bimaximal [5] 
matrix when e — rj — (5 = ^ = 0. If we allow the param- 
eters (e,7;,/3,^) to take on arbitrary values, the matrix 
above can be used as a general parametrization of the 
mixing matrix. Also, it is related to the familiar "stan- 
dard parametrization" by — e^, 3^2 — ^ + 
S'lg = 5 + 5Trj2) and /3 is a complicated function but 
/3 ~ ^C^5i3 for (e,?7,0 < 1- From VFo, we find read- 



1/6, 0:30 ~ 0, and Xio + ym 



ily xio ^ 1/3, X20 
(* = 1,2,3). 

In the flavor basis, the effective Hamiltonian for neu- 
trinos propagating in matter is given by -ffgff 
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(11) 



where rni, m2, and 7713 are the neutrino masses in vac- 
uum, Vq is the mixing matrix in vacuum, and the induced 
mass A = V2GFneE. 

The matrix H can be diagonalized, 

[D, \ 
H = VDV"^ = V\ D2 \V\ (12) 

where Di = Mf is the squared mass in matter. To study 
how the elements of V evolve in matter, one may start 
with dH/dA, 



dH _ d 
11' dl 
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[VDV^] = 



(13) 



unknown element [{dV^ / dA)V]ii is rephasing dependent, 
and can be set to vanish. 

While Eq. (|16p is rephasing dependent, it may be used 
to compute rephasing invariant quantities, e.g., 
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.dJ 
''dA' 



(17) 



After some algebra, separating the real and imaginary 
parts, in addition to using different T^^^s, we obtain the 
evolution equations for all {xi,yi) and d\nJ/dA, which 
are collected in Table I. The evolution equations obtained 
here are entirely analogous to the familiar RGE of mass 
matrices. In both cases, the effective Hamiltonian con- 
tains a parameter, the energy for RGE, and A for neu- 
trino propagation. The respective evolution equations 
can be used to solve for eigenvalues and mixings as func- 
tions either of the energy scale, or of A. 

The symmetric form of these equations allows us to 
find readily the result: 



D2){D2-D3)iD3-Di)]=0, (18) 



i.e., the product [J{Di - D2){D2 - D^){D:i - Di)] is a 
constant as A chaiiges, a well-known result derived with 
different methods [9|. 

In addition, by writing down the evolution equations 
for [d/dA) ln(xi — 2/i), from Table I, we find another "mat- 
ter invariant" : 



J2 



dA [xi - yi){x2 - y2){x3 - ^3) 



0. 



(19) 



Or, [JV(|Vii|2|\/i2|2|yi3|2)] = constant. When we 
use the "standard parametrization" , it is seen that 
JVi\Vii\^\Vi2\^\Vi3\^) = SplA ^-e-, 5j^in2023 is in- 
dependent of A, a result obtained earlier [ll]. 

The evolution equations for (a;, y) also have a structure 
akin to that of the fixed point of single variable equations. 
It can be verified that, if Xi + yi = {i = 1, 2, 3), then 



Eq. is then sandwiched by V'' and V, 



dA dA dA 



\Vn\^ Vi2V,\ VuV,\ 

ViiV,*2 1^12 P ^13^1*2 



ViiV.l V12V, 



13 



\Vl3\' 



(14) 

Taking the diagonal and off-diagonal terms of Eq. (ITil) . 
and following the procedures in Ref.f7|, we find 



— = \Vu\ ^Xi- yi. 



(* = 1,2,3) (15) 



dV, 



dA 
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Eq. m follows from [{dVydA)VU = V*V,k/{D,-Dk), 
i ^ k, which can be inverted to solve for dV/dA since the 
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0, 



J = (1,2,3). 



(20) 



This result is understandable since the conditions Xi + 
yi — are equivalent to W2i = W^i , which, in turn, imply 
that the effective Hamiltonian H has a fj, — t exchange 
symmetry [sj . This symmetry is clearly independent of A 
in Eq. ([TlT) . resulting in Eq. ([20)1 . Note also that there are 
actually only two independent constraints in Xi + yi = 0. 
Given any two of them, say for i = 1,2, we can use Eq.(4) 
to derive x^+ys = 0. Thus, the set of evolution equations 
has a "fixed surface", points on the surface defined by 
Xi + yi — stay on it as A varies. 

While analytical solutions to the equations in Table 1 
are not available, as we will see, given the known physical 
parameters, one can exploit certain characteristic prop- 
erties thereof to arrive at simple, but fairly accurate, so- 
lutions to these equations. 
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l/{Di ~ Da) 


1/(D2 - D3) 


1/(L»3 -Di) 


dx\/dA 


X1X2 - 2x17/2 + yiJ/2 


-X1X2 + X1X3 + 2/12/2 — yiy-j, 


-X1X3 + 2xi2/3 - 2/12/3 


dx2 /dA 


-X1X2 + 2x2yi - yij/2 


X2X3 - 2X22/3 + 2/22/3 


X1X2 - X2X3 - 2/12/2 + 2/22/3 


dx^/dA 


-X1X2. + 2:2X3 + J/iJ/s - 1/21/3 


— X2X3 + 2X32/2 " 2/22/3 


X1X3 - 2x32/1 + 2/12/3 


dyi/dA 


-X1X2 + 2x2yi ~ yiy2 


-X1X2 + X1X3 + 2/1J/2 — 2/12/3 


X1X3 - 2x32/1 + yi2/3 


dy2/dA 


X\X2 ~ 2xiy2 + J/lJ/2 


--X2X3 + 2X32/2 — 2/22/3 


X1X2 — X2X3 — 2/12/2 + 2/22/3 


dy-i/dA 


-X1X3 + X2X3 + yij/3 - 2/22/3 


X2X3 - 2X22/3 + 2/22/3 


-X1X3 + 2xi2/3 - 2/12/3 


d{\nJ)/dA 


— Xi + X2 + J/l — 2/2 


-S2 + X3 + 2/2 — 2/3 


Xl - X3 - 2/1 + 2/3 



TABLE I: dxi/dA, dyi/dA, and d(ln J)/dA are expressed as sums of terms in l/{Di — D2), 1/(^2 — D3), and l/{Ds — -Di). 



Experimentally, it is known that Sq = — m\ = region, the explicit, approximate, solutions are 



7 X lO-^ey^ Ao 



3 X 10"3ey2 



that (5o/Ao <C 1 (we assume the "normal" ordering of 
neutrino masses. The "inverted" case can be similarly 
treated). Note that these values are relevant to long 

baseline experiments since A = y/2GFneE ~ (7.6 x 
10"5eT/2)(^/GeF)(p/gcm-3). 

Since 5q <C Aq, we expect that the three- flavor prob- 
lem can be approximated by a pair of well separated 
two- flavor problems Indeed, the structure of the 

differential equations in Table I shows that the vari- 
ables (xi,yi) evolve slowly as a function of A except 
for two regions, where Di « D2 and D2 ~ -D3, cor- 
responding to the two resonance regions. More pre- 
cisely, let us denote by {Aq, Ai, Ai, Ah, Ad) the values 
of A in vacuum {Aq ~ 0), at the lower resonance 
{Ai,[d{Di — D2)/dA]Ai = 0), intermediate range {Ai), 
higher resonance {Ah, [d{D2 — D3)/dA]A,^ — 0), and for 
dense medium {Ad). Rapid evolution for {xi,yi) only 
occurs for Ak Ai and A^i Ah- 

For < A < Ai, we need only to keep the terms cx 
l/{Di- D2) in Table I. It is seen that 



d{xi - y2) _ d{x2 - Vi) _ d{x3 - 2/3) 



dA 



dA 



dA 



0. (21) 



S = [2biA^-AciA + Sl]'^^, 
Xl = ^[{xio + X20) - {bi A- ci)/S], 

X2 = ^[{xiQ + X20) + {bi A -ci)/d], (24) 

with ci = So{xio - X2o)- Also, xi + yi ^ X2 + y2 = 0, 
= 2/3 = 0. From Wq, we have 6/ = 2(a;io-l-a^2o)^ — 1/2- 
We see thus, as A sweeps through the lower resonance 
region, S goes through a minimum, Xi decreases and X2 
rises while keeping xi + X2 — 1/2. After the resonance, 
for large A (A > 60), S ~ A, xi ^ 0, and X2 -> 1/2. 

A similar analysis can be done for the region Ai < 
A < Ad- Here, the starting values {Ai <^ A < Ah) are 
Xl ~ 2/1 ~ 0, a;3 ~ 2/3 ~ 0, X2 1/2, 2/2 ^ -1/2. 
The differential equations are dominated by terms pro- 
portional to 1/{D2 — -D3), and they satisfy 



d{xi - 2/1) _ d{x2 - 2/3) _ d{x3 - 2/2) 



dA 



dA 



dA 



= 0. 



(25) 



With A = D3 — D2, the approximate equations near 
A^ Ah are then 



Given Wq (Eq. ^), with a;« + yio = 0, Eq. ^ yields 
Xi + yi = 0. Thus, we expect that < A < Ai, 
while the individual variables {xi,X2,yi,y2) are rapidly 
changing, X3 ~ 2/3 — ^i^) stay small, as do the combi- 
nations Xl + yi = X2 + y2 — 0, and xi + X2 — constant, 
yi + y2 — constant. The differential equations can then 
be approximated, with S = D2 — Di, by 

dxi ^ -4xiX2 ^ dx2 dS ^ 



It follows that 

_d_ 
dA 
d_ 
dA 



\x1X2S ] = 0, 



[{xi-X2)5]^~2{xi+X2)'' = -bi, (23) 



where bi ^ 2(a::io -t- X2o)'^- So, in the lower resonance 



dx2 ^ ~'^X2X3 ^ dxz dA 

together with 2:2 + J/2 — 2^3 + 2/3 — 0, while xi and 2/1 are 
slowly varying so that a;i ~ 2/1 — throughout. 
The solutions are 

A=[2bhA' -AchA + Al]'/', 

Xl = -^{x2a + 2:30) - {bhA - Ch)/A\, 

X3 = ^[(3^20 + 230) + {bhA - Ch)/A], (27) 

with bh = 2{x2o + 2:30)^ and Ch = Ao(x2o - a^so)- 

Thus, as A goes from Ai through Ah to Ad, the changes 

for {xi,yj) are: X2 ^ 1/2 0; 2/2 ^ -1/2 0; X3 ~ 

0^ 1/2; and 2/3^0^ -1/2. 

Our results can be summarized by the matrices W at 
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FIG. 1: Numerical (solid) and approximate (dot-dashed) solutions for (a) all Ds{A), D2{A), and Di{A), and (b) the enlarged 
plot of D2{A) and Di{A) in < A/ So < 10. 




FIG. 2: The numerical (solid) and approximate (dot-dashed) 
solutions for xi{A), X2{A), and X3{A). Note that yi{A) ~ 
-x,{A). 



A = {Ao,Al,A^,Ah,Ad)■. 




X3 

. 5t 
0.4 
. 3 
0.2 
. 1 



A/6o 



10^^ IQ-^ IQ-^ 1 10 10^ 10^ 



A/So 



''-^ 10"^ 10"^ 1 10 10^ 10^ 



, A/6o 



Wn = 



2/3 1/3 

1/6 1/3 1/2 

1/6 1/3 1/2 

/ 1 

1/2 1/2 

^ 1/2 1/2 

/ 1 

1/2 1/2 

y 1/2 1/2 



W, = 



/ 1/2 1/2 

1/4 1/4 1/2 
^ 1/4 1/4 1/2 

1/2 1/2 
1/2 1/4 1/4 
1/2 1/4 1/4 



(28) 



Together, these matrices exhibit the remarkable simplic- 
ity of the PMNS matrix as A varies from to oo. Note 
that all of the matrices have at least one zero, Wu ~ 0, 
implying xi — yi — 0. Also, they have equal elements in 
their second and third rows, W2i = W^i, so that Wu = 
or Xi + Hi = 0. As a consequence, using the unitarity 
conditions, the W matrix is completely fixed by its el- 
ements in the first row, Wu. These elements, in turn, 
control dDi/dA, Eq. (fT5|) . Thus, the progression of W as 
a function of A can be read off from the plot of Di{A), 
which is given in Fig. 1. 

It is straightforward to numerically integrate the evolu- 
tion equations for (x, y). To do this we choose the initial 
values (in Wq) e = 0.17, /? = 0.02, corresponding to the 
experimental bounds IKsP < 0.03 [1] and an assumed 
CP violation phase cose/) =1/4. Also, ^ = = 0. In Fig. 
2 the results are compared to the approximate solutions 
obtained earlier (Eqs. (1^ and (P7|). The agreements 
are quite good. The evolution of is shown in Fig. 3. 
Compared to its vacuum value, it is seen that, except for 
some enhancement near A = Ai, tends to decrease 
with increasing A^ as one would expect from Eq. (jl8|) . 

In conclusion, in this work we derived the evolution 
equations for the neutrino parameters as a function of 
matter density. We found analytic, approximate, but 
simple solutions of these equations for values centered 
around the known neutrino parameters. This is possible 
because of two fortuitous circumstances: 1) the neutrino 
mass differences are widely separated, enabling one to use 
the two-flavor resonance approximation; 2) the mixing in 
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FIG. 3: The evolution of from the numerical (solid) and 
the approximate (dot-dashed) solutions. 




A/5o 



vacuum satisfies lio+Z/io — 0; which happens to he on the 
"fixed surface" of the evolution equations, so that Xi + 
Hi ^ for all A values. These solutions are summarized 
in Eq. (|28p . vifhich exhibits the extraordinary simplicity 
of 14^ as a function of A. These results are found to be 
quite accurate when we compare them to those obtained 
by numerical integration of the equations. It is hoped 
that they will be useful in assessing the matter effects in 
connection with the long baseline experiments, as well as 
efforts to determine CP-violation in the leptonic sector. 
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